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The precession of the HH 111 flow in the infrared 

Noriega- Crespo, A.^, Raga, A. C. ^, Lora, V.'^, Stapelfeldt, K. R.'* and Carey, S. J.^ 

ABSTRACT 

We present Spitzer IRAC images of the HH 111 outflow, that show a wealth 
of condensations/knots in both jet and counterjet. Studying the positional dis- 
tribution of these knots, we find very suggestive evidence of a mirror symmetric 
pattern in the jet /counterjet flow. We model this pattern as the result of an 
orbital motion of the jet source around a binary companion. From a fit of an 
analytic, ballistic model to the observed path of the HH 111 system, we find that 
the motion in a binary with two ~ 1 Mq stars (one of them being the HH 111 
source), in a circular orbit with a separation of ~ 186 AU would produce the 
mirror symmetric pattern seen in the outflow. 

Subject headings: circumstellar matter — stars: formation — ISM: jets and 
outflows — infrared: ISM — Herbig-Haro objects — ISM: individual objects 
(HH 111) 

1. Introduction 

Over the past two decades considerable progress has been made in understanding the 
formation of stellar jets and their influence on the surrounding interstellar medium (see, e. g.. 
Bally, Reipurth & Davis 2007; Arce et al. 2007). Nevertheless, some fundamental questions 
remain unanswered, mostly related with the nature of the acceleration and collimation of 
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the jets themselves, the properties of the 'launch region' and how both of these processes can 
be influenced by nearby objects, since many outflows arise from binary or multiple young 
stellar systems. 

Since the discovery of HH 111 (Reipurth 1989) using standard narrow band imaging 
techniques to map the outflow in Ha and [S II] AA 6717/6731, HH 111 has become an icon of 
highly coUimated stellar outflows, and has been the subject of a number of studies across the 
entire wavelength spectrum to determine the proper motions of their 'knots' (Hartigan et al. 
2001; Coppin, Davis & Micono 1998), kinematics (Reipurth, Raga & Hcathcotc 1992; Raga 
et al. 2002), ionization structure (Reipurth 1989; Davis, Hodapp & Dcsochcs 2001; Nisini 
et al. 2002), circumstcUar environment (Stapelfeldt & Scovillc 1993; Yang et al. 2007), 
shock properties (Noriega-Crespo, Garnavich & Raga 1993; Morse et al. 1993), molecular 
flow morphology (Cernicharo & Reipurth 1996; Lefloch et al. 2007), and ejection history 
(Masciadri et al. 2002; Raga et al. 2002), among others. The original optical images were of 
the West lobe of HH 111, and the Eastern counterjet was later identifled in the near infrared 
(NIR) (Gredel & Reipurth 1994; Davis, Mundt & Eisloffel 1994). 

HH 111 hes in the L1617 cloud in Orion at a distance of ~ 417pc and is oriented almost 
E-W (PA~97.5°) with a length of ~ 11' in our IRAC images (see Fig 1). HH 111 is part 
of a quadrupolar flow (Reipurth ct al. 1999; Rodriguez et al. 2008), with HH 121 arising 
nearly perpendicular (PA~27.3°) to the HH 111 flow, suggesting that the VLAl source is 
part of a close binary system. Reipurth et al. (1999) proposed that this binary could have 
a separation of ~ 0".l (~ 50 AU). Rodriguez ct al. (2008) detect two elongated structures 
with a separation of ~ 0".036 (~ 15 AU), which could correspond to the previously proposed 
binary. 

Finally, we should note that HH 111 is the central region of a "giant HH jet", with two 
lobes extending out to ~ 3.5 pc from the outflow source (Reipurth et al. 1997). The far 
ends of these lobes (HH 311 and HH 113, sec Reipurth et al. 1997) show a point-symmetric 
deviation from the direction of the inner HH 111 outflow of ~ 4°. This change in direction 
between the inner and outer region of the HH 111 outflow is evidence for a precession of 
the outflow axis, with a period Tp 3.5 pc/150 km s~^ ^ 23000 yr (considering the proper 
motion of ~ 150 km s~^ measured by Reipurth et al. 1997 for HH 311 and 113). 

Recently, we have shown (Raga et al. 2011) that by using observations from the Spitzer 
Space Telescope (Spitzer) (Werner et al. 2004) in the mid-infrared, where the extinction 
is considerably lower than in the optical and near infrared, one can use the jet /counter jet 

symmetry to set some of the strongest constraints yet on the size of the jet formation region. 
In this study, we take advantage once again of the superb sensitivity and angular resolution 
of the Infrared Array Camera (IRAC) (Fazio et al. 2004) on board of Spitzer to study the 
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positional distribution of the "knots" along the Herbig-Haro (HH) 111 flow, finding mirror 
symmetric deviations from the jet/counterjet axis. This is the first time that such a mirror 
symmetric pattern is seen in an HH outflow, which can be interpreted in terms of the presence 
of a binary outflow source. 

2. Data 

The observations of HH 111 are part of our original Spitzer General Observer (GO) 
program 3315 (PI Noriega-Crespo) obtained with IRAC and the infrared photometer MIPS 
(Rieke et al. 2004) in March 28, 2005. The data have been downloaded from the Spitzer 
Legacy Archive and the quality of the final images (Post Basic Calibrated Data or Post- 
BCD; S18.7 products) is outstanding, so that no further processing was required. The 
HH 111 outflow is part of a larger mosaic of the L1640 region obtained with IRAC in its four 
channels (1, 2, 3, 4) = (3.6, 4.5, 5.8 & 8.0 /xm) covering a FOV of ~ 19 x 23' (the result 
of a 5x4 array map with a 260" stepsize) and with a total integration time per pixel of 360 
sec, using 30 sec High Dynamic Range (HDR) exposures. The final images are sampled with 
0.6" per pixel, nearly one third of the standard ~ 2" IRAC angular resolution. 

The IRAC images (Fig 1.), in comparison with the published NIR images (e.g. Davis, 
Mundt & Eisloffel 1994; Fig 5), display a similar morphology, but given their higher signal- 
to-noise and lower extinction permit us to define better the knots and identify a couple more 

on each side close to the VLAl source. Like in other protostcllar outflows observed with 
Spitzer (e.g., Noriega-Crespo et al. 2004; Looncy, Tobin & Kwon 2007; Tobin et al. 2007; 
Ybarra & Lada 2009; Raga et al. 2011) HH 111 is brighter in channel 2, although it is well 
detected in all 4 IRAC channels. The emission is most likely dominated by the H2 pure 
rotational lines (De Buizer & Vacca 2010). 

In the following we describe the emission peaks as condensations and/or knots. It is clear 
from the HST images at high angular resolution that many of these structures arc bowshock- 
like (working surfaces) and that they can overlap one another. We use the intensity peaks 
to define the HH 111 fiow. The measurements of the positions of the condensations/knots 
were carried out using the IRAF imexam routine interactively on the IRAC channel 2 image. 
For compact "knots" a two-dimensional gaussian was used for determining their centroids; 
for the more diffuse structures, we use the position of the pixel with the peak emission. The 
uncertainty of the centroids ranges from 0.1 to 0.3 of a fraction of a pixel for the bright and 
fainter knots, respectively. Since the "knots" are compact in most cases, we measured the 
fiux densities using a circular aperture and an aperture correction. For the compact "knots" 
we used a 2.4" radius for the photometric measurements with a 1.221 aperture correction; the 
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background was measured off the jet/counterjet ffow axis. For 7 "knots" that were obviously 
more extended or irregular than this aperture (6.23 x 10^^° steradian comparable 
polygonal area was used without an aperture correction. Some of these "knots" do have a 
low surface brightness, but it was still possible to measure the flux density and centroid. 

In Table 1 we give the relative positions in X (along the outflow axis) and Y (perpen- 
dicular to the outflow axis) with respect to the central source, VLAl at RA=05:51:46.25 and 
Dec=02:48:29.5 (Reipurth et al. 1999) of the "knots" seen along the jet and the counterjet. 
Table 1 contains also their RA and Dec coordinates, flux density and uncertainty (both in 
mJy). 



Figure 2 shows the inner region of the HH 111 outflow (top) with the measured (X, Y) 
positions of the knots with respect to the VLAl source. The bottom plot shows the knot 
positions with an expanded F-coordinate. It is clear that both the jet and the counterjet 
show an excursion in the — y-direction at a distance of ~ 40" from the source, and an 
excursion in the +y-direction at ~ 115" from the source. These excursions could be part of 
a mirror symmetric spiral shape with a step of ~ 2 x (115" — 40") = 150". 

These mirror-symmetric excursions of the jet/counterjet system could possibly be the 
effect of an orbital motion of the outflow source. Ballistic flows from sources in circular and 
elhptical orbits were described by Masciadri & Raga (2002) and Gonzalez & Raga (2004), 
respectively. 

For a ballistic flow, the path of the beam of a constant velocity jet ejected parallel to 
the orbital axis (of an outflow source in a circular orbit) is given by : 



where x is the axial coordinate, {y, z) are axes parallel to the orbital plane, ip is the orbital 
phase. To is the orbital period and k = Vo/vj is the ratio between the orbital and the jet 
velocities. As can be seen from a comparison of Equations ([I]l2]) with the results presented 
by Masciadri & Raga (2002), in these equations we have considered that the orbital radius 
is negligibly small compared to y and z. 



3. The mirror symmetry of the HH 111 jet/counterjet system 




(1) 




(2) 
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From Equations ([TH2]), we see that the spiral described by the jet beam has a step 

L = VjTo . (3) 

From Equation ([1]), we calculate the predicted locus of the jet /counter] et structure, and 
project it onto the plane of the sky, assuming an angle of 10° between the orbital axis and 
the plane of the sky (Reipurth et al. 1992), and that the y-axis is parallel to the plane of 
the sky). We then carry out a least squares fit in order to obtain the values of the model 
parameters that best fit the HH 111 jet/counterjet locus within 150" from the outflow source. 
The results of the fit are shown in the left panel of Figure 2. 

The following parameters are obtained from the fit : 

1. an orbital phase = (101 ±38)°. This result implies that the line joining the two stars 
of the binary outflow source currently lies close to the line of sight, 

2. a spiral step L = (216 ±39)". For a distance of 414 pc to HH 111 this corresponds to a 
distance of (1.34±0.24) x 10^^ cm. Also, from the proper motions measured by Hartigan 
et al. (2001), we see that the knots at ~ 40" along the jet (corresponding to the first 
sideways excursion in the jet path, see Figure 2) have a spatial velocity Vj ~ 240 km s^^. 
From Equation we then obtain an orbital period = 1801 ± 317 yr, 

3. a K = 0.013 ± 0.008 (see equations [T]12]) . For a jet velocity of 240 km s~^ (see above), 
this implies an orbital velocity Vo = 3.1 ± 1.9 km s~^. 

We now use the relation (valid for a circular orbit) between the mass of the jet source 
Ml, the mass of the companion M2 = aMi, the orbital period To and the orbital radius (of 
the jet source) : 

' - ° " - (0.10, 2.02, 6.97) M0, (4) 



(l±a)2 2ttG 

where the three values in the third term correspond to (k, L) = (0.005, 177") (the lower 
limits for k and L, see above), (k, L) = (0.013, 216") (the values that result in the lowest S^) 
and {k,L) = (0.021,255") (the upper limits). 

Each of the three values of the third term of Equation corresponds to a curve in 
the M2 vs. Ml plane, which are shown in Figure 3. Together with the M2 < Mi condition 
(i.e. assuming that the source of the strong HH 111 outflow is the more massive star of the 
binary), we then obtain an allowed region for the masses of the HH 111 source binary. 

From Figure 3, we see that the observed shape of the HH 111 jet/counterjet system is 
consistent with a binary with two ~ IM© stars. Such a mass is consistent with the masses 
of T Tauri stars. 
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If we consider Mi = M2 = M pa 1 (consistent with the curved shape of the HH 
111 system, see Figures 2 and 3) and = 1800 yr, from the equations of a circular orbit we 
then obtain a separation 2ro ^ 186 AU between the two stars in the binary. Given that the 
orbial phase derived from the fit imphes that the stars in the binary are closely aligned with 
the line of sight, the projected separation between the stars would be substantially smaller 
than 186 AU. 



4. Summciry and conclusions 

We have presented new IRAC images of the HH 111 outflow. These images show 
chains of aligned knots along the jet and the counterjet. We have carried out astromctric 
measurements of the positions of these knots, and have also determined their fluxes in the 
IRAC 2 (4.5 fim) channel. 

The jet /counterjet knots within ~ 2' from the outflow source show mirror symmetric 
deviations from the outflow axis, which can be interpreted as the effect of an orbital motion 
of the outflow source. If we flt the observed shape of the projected path of HH 111 with a 
ballistic jet model (Masciadri & Raga 2002), we flnd that the observations imply the existence 
of a binary with a Tq = (1801±317) yr period, and an orbital velocity Vg — (3.1±1.9) km s~^. 
This range of periods and orbital velocities would be consistent with a binary formed by two 
~ 1 M0 stars, with a separation of ~ 186 AU. The model fit also imphes that the stars are 
currently quite closely aligned with the line of the sight, so that this separation might be 
consistent with the considerably lower projected separation implied by the observations of 
Rodriguez et al. (2008). 

Interestingly, the orbital period that we determine (tq ~ 1800 yr, see above) is ~ 13 
times smaller than the precession period Tp implied by the point-symmetrical "HH 111 giant 
jet" (see Reipurth ct al. 1997 and the discussion at the end of section 1). The resulting 
Tp/T"o ~ 13 ratio is consistent with the predictions from the precessing accretion disk models 
of Terquem et al. (1999). 

Raga et al. (2009) pointed out that the models of Terquem et al. (1999) implied that 
outfiows from binary sources would have mirror symmetries close to the outfiow source (due 
to the orbital motion) and point symmetries at larger distances (if the ejection direction 
follows the disk precession). Our new observations show that the HH 111 jet /counterjet 
system has mirror symmetry within ~ 2' from the source. This result, combined with the 
point symmetry of the HH 111 giant jet (Reipurth et al. 1997) is the first evidence ever of 
an astrophysical jet with the predicted fiip in types of symmetry between regions close to 
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and far away from the outflow source. 

It is evident that the mirror symmetric, side-to-side excursions of the HH 111 jet /counter jet 
are a quite subtle effect. If our interpretation of these deviations were incorrect, we would 
not expect that a comparison with a model of an outflow from an orbiting source would give 
us reasonable masses for the binary components. However, we find that a binary with two 

~ IMq stars and a separation of ~ 200 AU produces the observed spiral shape of HH 111. 
This completely reasonable result is a strong indication that we might indeed be interpreting 
the data correctly (see also the case for HH 211 by Lee et al. 2010). 

This work is based in part on observations made with the Spitzer Space Telescope which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA 
contract 1407. The work of AR and VL was supported by the CONACyT grants 61547, 
101356 and 101975. 
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Fig. 1. — The HH 111 flow in the mid-IR using the 3.6 (blue), 4.5 (green) and 8/im (red) IRAC 
channels. Notice the scattered hght (blue) near the source and the nearly perpendicular 
bipolar outflow HH 121. 
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Fig. 2. — The right panel shows the central region of HH 111 at 4.5/xm rotated so that its 
flow axis is vertical, and with the positions of the knots are shown with the white circles. 
The left panel shows the same knots positions (shown with crosses) with an expanded Y- 
axis (perpendicular to the outflow direction). The solid curve corresponds to a /t = 0.013 
ballistic jet (best flt model), and the dashed lines correspond to k = 0.005 and 0.021 (the 
three models have L = 216" and (p = 101°, see the text). The jet path predicted from the 
best fit model is also shown on the right hand side plot (white curve). 
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Fig. 3. — Parameter space for the primary and secondary masses of the binary system (Mi 
and M2, in units of Mq) derived from the curved shape of the HH 111 jet /counter] et system. 
The two curves represent the hmits for the allowed masses derived from fitting the ballistic, 
orbiting source jet model to the observed trajectory (see equation H]). The shaded region 
shows the parameter space allowed by the fit. The solid curve is deduced from the parameters 
of the best fit. If one assumes that M2 < Mi (i. e., that the source of the HH 111 outflow 
is the more massive star in the binary), the masses of the binary should lie beneath the 
M2 = Ml line (the straight diagonal cutting across the plot). 
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Table 1. Counter Jet (East) & Jet (West) Positions & Flux Densities 
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Table 1 — Continued 
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^Over a 1.058 x 10~^° steradian area 
^Over a 3.976 x 10~^° steradian area 
"^Ovcr a 6.431 x 10~^° steradian area 
'^Over a 8.292 x 10"^° steradian area 
^Over a 6.685 x 10~^° steradian area 
^Over a 7.615 x 10^^° steradian area 
^Over a 8.377 x 10~^° steradian area 



